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Abstract 

We study a model of competition for resource through a chemostat-type model where species 
consume the common resource that is constantly supplied. We assume that the species and resources 
are characterized by a continuous trait. As already proved, this model, although more complicated 
than the usual Lotka-Volterra direct competition model, describes competitive interactions leading 
to concentrated distributions of species in continuous trait space. 

Here we assume a very fast dynamics for the supply of the resource and a fast dynamics for 
death and uptake rates. In this regime we show that factors that are independent of the resource 
competition become as important as the competition efficiency and that the direct competition 
model is a good approximation of the chemostat. 

This is performed through asymptotic analysis, introducing different scales for the resource 
renewal rate and the uptake rate. The mathematical difficulty relies in a possible initial layer for 
the resource dynamics. 

The chemostat model comes with a global convex Lyapunov functional. We show that the 
particular form of the competition kernel derived from the uptake kernel, satisfies a positivity 
property which is known to be necessary for the direct competition model to enjoy the related 
Lyapunov functional. 

Key- Words: Ecological competition for resource; Direct competition; Multiscale analysis; Lyapunov 
functional, Adaptive dynamics. 

Mathematics Subject Classification: 34E13, 35K57, 47G20, 92D15, 92D40 



1 Introduction 

Evolutionary origin of species diversity has been one of the most important problems in evolutionary 
ecology. Individuals with similar traits (e.g. body size or shape) face strong competition, which might 
lead to the adaptive radiation. As a model of such competition-driven speciation, several theoretical 
studies have shown that the species distribution in continuous trait space will often evolve toward single 
peak or multiple peaks that are distinct from each other (see [181 QUI El EJ [IS] ) . Although many 
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previous results are derived from a simple model that assumes direct species competition of Lotka- 
Volterra type, competitive interaction among species generally occurs in competition for resource such 
as prey or nutrient. For example, birds with similar beak shapes are in competition because they 
utilize similar food resource. Another example is found in ecological stoichiometry where consumer 
species with similar C:N:P (carbon: nitrogen: phosphorus ratio) requirements experience competition 
for nutrient supply with their optimal C:N:P ratio (see pj5]). Thus, the competitive interaction should 
be modeled not directly but implicitly through the resource. 

Recently, we have proposed such a mechanistic model of species competition where species (with 
continuous trait) compete for common resource (which is characterized by another continuous trait) 
which is constantly supplied (chemostat-type model). We have shown that evolutionary branching, 
which is parallel to adaptive radiation, is possible in this chemostat-type model as well, supporting 
the results previously obtained by direct competition models. 

In order to faithfully capture the dynamics of species competition, a detailed mechanism such as 
species-resource dynamics should be explicitly modeled. On the other hand, specifying a particular 
mechanism narrows the applicability of a model. This is why conceptual and simple Lotka-Volterra 
type models of direct competition are widely and intensively studied. Then a natural question arises: 
Can we reduce a species-resource dynamics model into a Lotka-Volterra model? If so, how are pa- 
rameters (e.g., intrinsic growth rates, competition coefficients) determined by those of the original 
mechanistic model? The last question is interesting particularly when we think of a direct competi- 
tion model in continuous trait space. If we naively assume a Gaussian function as competition kernel, 
then it promotes the coexistence of infinitely many species (continuous distribution over continuous 
trait axis). But as soon as the kernel deviates from Gaussian, the evolution toward single or multiple 
peaks occurs, which is parallel to the emergence of distinct species [12]. The importance of this estab- 
lished result would increase if the conversion of resource-consumption kernel into direct-competition 
kernel does not guarantee that a Gaussian function is converted into another Gaussian function. 

This paper aims the reduction of a mechanistic model of species-resource dynamics to more concep- 
tual model of direct species competition. In Section [21 we describe our chemostat model of species- 
resource dynamics and a brief summary of our reduction result. Section [3] is devoted to a rigorous 
description of this result and the proof. In section U we extend our result to include a system with 
mutation. In section [U we study the Lyapunov functional property which is useful for evolutionary 
stability analysis. In section [6l we show some numerical examples to study how the reduced dynamics 
approximates the original dynamics. 



2 Model and a brief summary of the results 

In general, creatures require several different resource to survive and reproduce. Some resource might 
be in a strong competition with others, some resource might be sufficiently supplied, and some re- 
source might not be in competition at all. In case of phytoplankton species, resources are light and 
nutrients. Although tremendous amount of light is supplied from the sun, the competition for light 
still exists. Nutrient supplied from the bottom is also a limiting resource, and this dynamics might 
have much slower timescale compared to the competition for light (but slower dynamics do not imply 
less importance). There might exist some features that determine growth or death rates independently 
of resource competition. 

Inspired by such phytoplankton ecology, our mechanistic model for a chemostat with continuous 
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resource supply and continuous consumer population trait is given by 
§in{x, t) = n(x, t) [b(x) - d(x) + / K(x, y)R(y, t)dy] , 

(1) 

^ §- t R(y, t) = m(y) [R in (y) - R(y, t)] - R(y, t) J K(x, y)n(x, t)dx, 

where n(x, t) > is consumer species density with trait x € R and R(y, t) > is the common resource 
density with trait y 6 R. We assume that a death rate d(x) depends on species trait x. Reproduction 
rate consists of two parts: b(x) denotes a basic reproduction rate independent of the competition for 
the common resource, and J K(x,y)R(y,t)dy denotes a reproduction rate coming from consuming 
the common resource R. The resource R is supplied into the system with a rate m(y) so that R(y,t) 
approaches Rm(y) > 0. These drive the system (a chemostat type model). Consumption for resource is 
denoted by a trait dependent uptake kernel K(x, y) > that defines how species with trait x depends 
on resource with trait y. 

All coefficients (b(x),d(x), K(x,y),m(y), i?j n (y); although they are functions of traits, we call them 
'coefficients') can play an important role in the full system. Since we are interested in the reduction of 
this model to a direct-competition model, we hereafter consider some limiting situation. We assume 
that the dynamics has the following three different time scales. 

• Very fast dynamics for the common resource supply. This imposes that, in the absence of con- 
sumer species, resource distribution R approaches R{ n . 

• Fast dynamics for resource-consumption and the resulting reproduction and death. This implies that 
competition for the common resource is a primary factor for the growth of consumer species popula- 
tion. 

• Slow dynamics. Population dynamics of n is also affected by a basic reproduction rate b(x) and a 
death rate d s \ aw (x), both of which are independent of the common resource. 

In order to include these assumptions in the equations, we scale accordingly the coefficients 

• m is of order e~ 2 , 

• K and d are of order e~ l , 

• d(x) = \ f K(x,y)R in (y,t)dy + d s i ow (x), 

• a(x) = b(x) - d s i ow (x). 

A death rate must be sufficiently large to counterbalance fast growth based on resource-consumption 
(otherwise, the population will grow infinitely). For this purpose, in the third equation, we divide a 
death rate into two parts. The first term represents fast dynamics for death events that compensates 
fast resource-consumption. The second term, d s \ ow (x), represents the remaining slow dynamics for 
death events independent of the competition for resource. 
These result in a new formulation of ([T]) that is written as 

§- t n £ {x,t) = n £ (x,t)[a(x) + \ J K(x,y)(R £ (y,t) - R- m )dy], 

(2) 

^R e (y, t) = ^ [R in (y) - R £ (y, t)] - \R £ {y, t) J K(x, y)n £ (x, t)dx, 
where new coefficients satisfy m = 0(1), K = 0(1), and a = 0(1). 
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We complement the system with initial data 

n £ (0, x) = n°(x) > 0, R £ (0, y) = R°(y) > 0, 

As the very fast and fast dynamics become infinitely fast, R should approach i2j n , but it is not 
clear which distribution n should approach. When R is not close to -Ri n , the reproduction based on 
resource-consumption is dominant (fast dynamics with order e~ 1 ) so the dynamics of n is primarily 
determined by R — i?; n . However, as R — R{ n approaches to zero, resource-dependent reproduction 
might become no longer dominant. 

In this paper, we will show that an approximation of this dynamics is given by the direct competition 
model 

d f 

—n(x,t)=n(x,t)[a(x) — / c(x, x')n(x' , t)dx'~\ , (3) 

where the binary competition kernel c(x, x') is related to the coefficients in ([2]). We will show that it 
is given by 

c(x,x')= [ K(x,y)^^-K(x',y)dy. 
J rn{y) 

This equation can be seen as a conversion formula of kernels. The reduced competition-kernel c(x, x') 
is very sensitive to R m {y)/m{y) as well as to resource-consumption kernel K(x',y), but not to Ri n {y) 
or m(y) separately. Note that this kernel is of order one (recall that resource-consumption kernel is 
of order e _1 before rescaling) so the reduced direct competition model describes the slow dynamics of 
n(x, t). We will later show that the inclusion of mutation in consumer species' trait does not essentially 
change the result. 

One biologically interesting point of the equation ([3]) is that the growth rate of consumer species in 
the absence of competition for resource, which is described by a(x), still has considerable effect on the 
long-term dynamics of consumer distribution n. This is because as R becomes very close to i?; n , the 
dominant factor (fast dynamics) becomes weaker and a basic growth rate a(x) becomes as important 
as the competition for resource. In this limit, we recover a Lotka-Volterra type competition dynamics 
from a mechanistic model of a chemostat type. It is natural that we do not have the reduced dynamics 
of the resource, R, because the reduction is possible when the detailed dynamics of the resource can 
be neglected and when the effects of competition for the resource can be put into the new competition 
coefficients in a Lotka-Volterra type system (i.e. c(x,x')). 

Another question of our interest is whether a Gaussian uptake kernel K(x, y) is converted to another 
Gaussian kernel of direct competition c(x,x'). This is not necessarily true. Take 

K(x,y) = e~^-y)\ R . m{y) = e ~Pv\ m{y) = \ 
Then, choosing 7 such that 7(2a + (3) = a 2 so as to complete a square in the expression 

x '-j = e -[ax 2 +a(x') 2 -j(x+x') 2 ] J g-[(2a+/% 2 -2a(x+x')y+7(x+x') 2 ] dy 

_ — [ax 2 +a(x') 2 — y(x+x') 2 ] I it 
~ C Y 2a+/3 ' 

we find that c(x, x') is not gaussian because it is proportional to 

e -[ax 2 +a(x') 2 --y(x+x') 2 ] _^ e -(a--y)(x-x') 2 w h en /3 ^ 0. 

Mathematically rigorous notations and results will follow in the next three sections. Section [3] is 
devoted to a rigorous derivation of this result. It goes through asymptotic analysis and compactness 
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estimates that allow us to show that n £ converges to n and thus n is the population density. In 
other words our theory is not a perturbation analysis on n £ while it is on R £ ~ R\ n . Mutations are 
not included in the model equations (jTJ) — d3J) . In section HI we show how to include them and which 
new rescaling this leads to perform. Another remarkable property of the chemostat is the Lyapunov 
functional property already used in [16] . The asymptotic form of this Lyapunov functional is related 
to that for the direct competition model in section [5] and in particular we explain why, in this context, 
the corresponding operator is always nonnegative. 



3 Derivation of the direct competition 

In order to state our result we need some assumptions 

n°, R°, R m a'n L°°(R), I \-^-n°\ dx < C u (4) 

J ox 

a,M '■= \\ a \\oo < oo, m = minm > 0, m = maxm < oo, (5) 

K e BUC(R x R]0, K M := sup K(x, y) < oo, (6) 

I— K(x,y)\ <C 2 , | —a(x)\<C 3 . (7) 
ox ox 

With these in mind, we establish that the limiting problem for ([2]) is the direct competition model 

Proposition 3.1 As e tends to in (|2|), R e converges strongly to R[ n in L 1 ((0, T) x M) and n £ 
converges strongly in C(M + ; L 1 (M)) to n, the solution to (J3]) with the initial data n° and the competition 
kernel 

c(x,x') = j K(x,y)^^-K(x',y)dy. (8) 



It is more technical to improve the space-time convergence on R e . Indeed it leaves place for an 
initial layer which is needed because, in general, its limit R[ n is incompatible with the initial data R®. 

Proof. To begin with, we re-write the equation on R £ in (|2|) as 

R £ (y,t) - R in (y) = -e^j 2 ^- [ K(x,y)n £ (x,t)dx £ —^-R £ (y,t). 

m(y) J m{y) dt 

We inject this to compute the growth rate for n £ 

^ j K(x,y)(R £ (y,t) - R in )dy = - J K(x, y) ~7~^~ J K(x' ,y)n £ (x' ,t)dx'dy + I £ (x,t) 

= — J c(x, x')n £ (x', t)dx + I £ (x, t), 
as we use c(x, x') as defined in ([8]) and with 

I e (x,t) = -el j ^^R £ (y,t)dy + j K{x,y) ^ y) -^ M f K(x' ,y)n £ (x' ,t)dx'dy. 
x Linear space of Bounded Uniformly Continuous functions 
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Therefore the equation on n £ can also be written 



d 

—n £ (x,t) = n £ (x,t) a(x) - I c(x,x')n e (x' ,t)dx' + I e (x,t) . (9) 



Of 

Throwing away the terms in e, we find formally (|3j). 

To justify rigorously this limit, we need (i) estimates for compactness, (ii) explain how we can pass 
to the limit. 

(i) Estimates. First notice that we have 

< R E (y, t) < R M (y) ■= m&x{R°, R in ) e L 1 n L°°(M). 

Next, we define and write from the equation on R £ in PJ 

M e (t) := J n £ (x,t)dx, N £ (t) := J \R e (y,t) - R in (y)\dy, 



d m^^\\R M (.)\\ 

dt y ' e 

And from the equation on n £ in ([2]), we deduce 



-N E (t) + =N £ (t) < KmM e (t) . 



\^-M £ (t)\<a M M £ (t) + ^N £ (t). 
dt e 

It follows from these two differential inequalities that M £ (t), eN £ (t) and ^ J * N £ (s)ds are uniformly 
bounded in e on each time interval [0, T] by a constant denoted by Mt- One step further, from the 
bound on M £ and from the equation on N £ , one also concludes that 

M> < m e -^ + !Mit KMlfr , « e [o. n. do) 

e e 721 

(m,) Compactness in t. From this, and going back to the equation on n £ , we conclude that, still for 
t€ [0,T] 



/ 



dn e (x,t) 



01 



N(0) 



dx < a A /C(T) + K M *-^-e-B*l* + K M C(T). (11) 



This gives us that n e is uniformly equicontinuous in C((0, T); . 

(raj Compactness in x. Dividing equation ([2]) by n e and differentiating in x, we find 

d 2 d 1 f d 

-^lnn £ (x,t) = — a(x) + - J —K{x,y) (R £ (y,t) - R in {y)) dy. 

From ([7|) and (|10p we obtain that, for some constant D\ 

, d 2 



dxdt 

It follows that 



lnn e (x,t)| < D\. 



\-^lnn £ (x,t)\ < |^lnn°(x)|+Z>ii, 
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/ 



— n e (x,t)|cix < / n , r — nJx, t) dx + D\t I n £ (x,t)dx. 
ox J rv ) {x) 



Moreover, from the first equation of ([2|) and (|10p we obtain that, for some constant D 2 , 

d 

\—lnn £ (x,t)\ < D 2 , 

and thus: lnn e (x, t) < lnn°(x) + D 2 t, which gives 

n E (x,t) 

Using the above arguments and (jl]), we conclude that, for some constant D(T) and < t < T, 

J \-^-n £ (x,t)\dx < exp(D 2 t) J \^-n° \x)\ dx + D x t j n £ (x, t)dx < d exp(D 2 t) + D^tM e {t) < D(T). 



Using this and step (ii) we conclude from the Arzela-Ascoli Theorem that, after extracting a subse- 
quence, n £ converges in (^((OjT);^ 1 ! 



(iv) Passing to the limit. It remains to pass to the weak (distribution) limit in the equation on 
n £ written as ([9]). Because of the strong convergence of n £ and (llOp . the only difficulty is to pass to 
the limit in the term 

. .9 / K(x,y) . . , d r . . f K(x,y) _ . . , -, 
en £ {x,t)— J R £ (y,t)dy = e-[n £ (x,t) J Re(y,t)dy\ 

f K(x,y) _ . , , d , . 

- £ J ^ R ^ t)dy ot n ^- 

The first term converges weakly to (multiplying by a test function, after integration by parts all the 
terms are bounded and multiplied by e). For the second term, we use the estimate (lllj) to reach the 
conclusion that it converges in L 1 to 0, thanks again to the multiplying factor e. 

Since uniqueness of weak solutions for the direct competition model is proved in [8], we conclude 
that the full family n £ converges and not only subsequences. □ 

4 Mutations 

In the present context of population equations structured by a physiological trait, there are several 
possibilities to represent mutations that have been widely used [3j \10\ [8] . Integral oprators or diffusion 
operators can be derived from stochastic individual based models; see for instance (5J [6]. Mathemat- 
ically they have many similar properties, in particular asymptotic analysis can be carried out using 
similar methods in the regime of small mutations for a long time of observation leading to specia- 
tion phenomena [21 dj [TJ] although stochastic individual methods are also used for the same purpose 

Here, we restrict ourselves to representing the mutations by a diffusion operator with intensity fi 2 , a 
very small rate. According to the mathematical theory in [10} 117]. it is natural to rescale time so that 
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the new timescale unit is 1/fi, very long timescale compared to the original 'generation' timescale. 
This leads to the competition model with mutations 

Hl%n e (x,t) = n e (x,t)[a(x) + ~ J K(x,y){R £ {y,t) - R m )dy] + fi 2 An £ , 

(12) 



li&Reiy, t) = ^ [R m (y) - R £ (y, t)] - \R e {y, t) f K(x, y)n £ (x, t)d 



Because both the uptake rate e and the mutation rate /x can be considered as small, it is interesting 
to see how the solutions behave in the different regimes of smallness. 

We can formally follow the lines of the analysis in section [3] and find direct competition with 
mutations as the limit of infinitely fast uptake rate 

d f 

fj,—n(x,t) = n(x,t)\a(x)— / c(x, x')n(x', t)dx'] + fi 2 An. (13) 
ot J 

One can also consider the small mutation rate and, assuming the initial data is well prepared, use 
the Hopf-Cole unknowns 

Ufj, = //ln(n £ ), u°=^ln(n°). 

This leads now to study the equation on the limit u £ (x,t) of u^(x,t) and on the limit measure n £ as 
[i — > 0. Following (10} I17j. one finds the constrained Hamilton- Jacobi equation 

d_ u — i 1 f „,\(r> („. +\ T>. \A*. i \9ujl\2 



dt 



a(x) + lJK(x,y)(R £ (y,t)-R in )dy + \^\ 2 , msixu £ {x,t) = Vt > 0, 



~7T- 



[R in (y) - R E (y, t)] - \R £ {y, t) J K(x, y)n £ (x, t)dx = 0, 



and n £ is supported in the zeroes of u £ (usually points). 
This can be further analyzed in the regime e small and leads to 

—u = a(x)+ / c(x,x')n(x' ,t)dx' + \ — \ 2 , max?x(x,t) = Vt > 0, 
ot J OX x£R 

still with the information that the measure n(t) is supported by the points where u(t) vanishes. This 
is also the limit of (1131) in the regime of small \i. 



We can conclude that the analysis of small mutations, both for the chemostat model and the direct 
competiton model, is compatible with the multiscale analysis in section [3] that relates both models. 
This is the regime of interest that we choose later in section [6] for numerical illustration. In this 
regime, we recover also, at a populational level, the results of standard adpative dynamics [TTJ [S] and 
the references therein. 



5 Lyapunov functionals for the two competition models 

The large time behaviour of solutions to (|2|) has been studied in [16]. It is proved that steady states 
cannot be globally positive (n £ , R £ ) and they have to be either Dirac masses (speciation) or continuous 
with a support of n £ that is small enough (except when the solution goes extinct that might happen if 
renewal/death rates are too high). When continuous, the steady state globally attracts all trajectories 
if they satisfy a particular sign property characterized the so-called Evolutionary Stable Distribution 
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(ESD in short), a notion that we recall below. The proof relies on a Lyapunov functional which has 
also been proved to exist for the direct competition model in [13] under the condition of positivity of 
a certain operator. 

Our interest here is to understand the relation between the two Lyapunov functionals, for ((2|) and 
([3]), and to understand why the operator positivity condition is important for Direct Competition but 
does not appear in the Competition for Resources. 

We first recall that an ESD is defined as 

Definition 5.1 (Evolutionary Stable Distribution for ([3]), [13J ) A nonnegative bounded measure 
n is called an ESD for the direct competition equation (J3j) 



a(x) — J c(x, x')n(x')dx' = \/x € Supp n, (14) 
« W -/c(„W)<fa'<0 V« R \Supp B . (15) 



v')n(x')dx' 

Definition 5.2 (Evolutionary Stable Distribution for , |13L I16j ) A state is (n £ ,R £ ) is called 
an ESD for the competition for resources equation ([2]) if n £ is a nonnegative bounded measure and 

a(x) + - J K(x,y)[R £ (y) - R in ]dy = Vx£Suppn e , (16) 
a{x) + ^ J K{x,y)[R £ (y)-R m }dy<0 VxGM\Suppn £ , (17) 

^[R £ (y) ~ Rin(y)} = R £ {y) J K(x,y)n £ (x)dx Vy G R. (18) 

For an ESD and following |16j . we define the Lyapunov functional for the model of competition for 
resources as 

S cr (t) = - J n £ (x) In n £ (x,t)dx - J R £ (y)ln R £ (y,t)dy + J n £ (x,t)dx + J R £ (y,t)dy. 
The dissipation associated with this functional is given by the inequality 

*t(t)-.= D cr (t)= -I ^tt) (My) - Re(y,t)) 2 dy ^ 

+ j n £ (x,t) (a(x) + IJ K{x,y)[R £ (y) - R in (y)]dy) dx < 0. 



Because we expect that, according to the proof of Theorem 13.11 the expansion holds R £ (y,t) = 
Rin(y) + 0(e), the limiting Lyapunov functional is reduced (up to an additive constant) to 

Sdc(t) = — J n{x) lnn(x, t)dx + J n{x,t)dx. 

This is indeed the Lyapunov functional used in for direct competition but -^Sddt) is negative only 
when a positivity condition holds true. This condition appears naturally in our framework because 
we can compute the dissipation associated with 5^ c (t) 



9 



A direct computation of D^c using ([3j) is of course also possible but is less instructive than an expansion 
in e of D cr that uses simply 

Re(y) = Rin(y) - £ -^v I K(x, y )n e (x)dx. 

m{y) J 

In the limit e — > 0, we obtain 

AfcW = -J^tfrU K(x, y) {n(x, t) - n(x))] 2 dy 

(20) 

+ Jn(x,t) (a(x) - J K(x,y)^- J K(x',y)n(x')dx'dyj dx. 
Using the expression of c(x,x') in (|8|) yields the usual dissipation for the direct competition model 
Ddc{t) = — / f c(x, x')(n(x, t) — n(x)) (n(x', t) — n(x'))dx dx' 
+ f n(x,t) (a(x) — f c(x,x')n(x')dx'^ dx. 

The particular form of c(-, •) in (jSJ) makes that the integral with the quadratic term in n — n is 
always positive (we referred to that as operator positivity) and the definition of an ESS makes that 
the linear term also gives a nonpositive contribution. 

6 Numerical illustration 

How good is the approximation of competition through resources by direct competition? Beyond the 
pure theoretical statement in section [3l this can be illustrated thanks to some numerical simulations. 
We have performed such comparisons in the case with mutations given by (| 12 [) . We assume that a 
resource uptake kernel and a distribution of resource supply are Gaussian. A trait-dependent (and 
competition-independent) growth rate takes the maximum at x = 0. We have used the coefficients 
given as follows 



K (x,y) = 7^= ex P( 7T2 — )' R in(y) = 7== ex P(-7r-2-)' a(x) = l-x 2 , (21) 



(x-y)\ T> / \ Mm ,_J? 

/rr—^f\ n 2 /' ttm^y) — /—— ex Pv n 

0-A-V27T 2,a R <7inV27T 2(7;, 

and m{y) is taken as a constant. See section[2]for the converted kernel of direct competition. Through- 
out this section we have taken cifc = a- m = .5, fj, = 0.005 and we vary e, m and M- m . 

The initial data n° is a gaussian centered at —0.8, the mean variance is denoted by n and R is ini- 
tially chosen equal to i?i n . With all our choices of parameters, the solution is inititally monomorphic 
and undergoes a dimorphic branching. Then it stabilizes as depicted in Figure [TJ 

We first show how well approximated is the solution to f)12|) by that of (|13p . To do so, choose 
m = M- m = 1 and choose two values of the parameter e, namely e = .1 and e = .001. For e = .1 
equilibrium of resource R ~ Ri n is not achieved since we obtained the relative error max ^p^ 1 " = .18 

but the case e = .001 can be considered as the direct competition case since max 1 R m = .0018. The 
population density are depicted in Figure [2] which show that the dimorphic situation is well kept in 
both solutions even though the total population is under estimated in the direct competition model 
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Figure 1: With our choice of parameters, the dynamics undergoes a branching from monomorphic 
to a dimorphic population. This figure depicts the typical isolines of the solution to (fT2|) (horizontal 
axis is x, vertical axis is time). 




Figure 2: Comparisons between the model of competition for resources (dotted curve) and the direct 
competition model (continuous line). Left: the dimorphic population density as a function of the trait 
x for large times. Right: total population J n(x,t)dx as a function of time. 
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Figure 3: These figures illustrate that the direct competition kernel only depends upon on the ratio 
R in /m and thus we can expect that the solutions to the model of competition for resource also does 
so. (Left) we compare for e = 0.05, the population density solutions with m = M- m = 1 (dotted 
line) and with m = M m = 10 (continuous line). (Right) we compare for e = 0.05 the solutions with 
m = Mi n = 1 (dotted line) and with m = 1.5, M- m = 1 (continous line). When the ratio is preserved 
(even though coefficients differ by a factor 10) the solutions are very close. But a small variation of 
the coefficients, that also changes the ratio, induces a visible difference on the solutions. 

by a factor of approximately 25 percent. 

The competition kernel ([8]) in the direct competition approximation has a remarkable property that 
it only depends upon the ratio R\ n /m. Therefore we can expect that the solution to (fT2|) depends 
only midly on this ratio. This is illustrated in Figure [3l We use the value e = 0.01, and compare the 
solutions with two values of m and M{ n and varying the ratio. 
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